Sensitivity of arrest in mode-coupling glasses to low-g structure 
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We quantify, within mode coupling theory, how changes in the liquid structure affect that of the 
glass. Apart from the known sensitivity to the structure factor 5(g) at wavevectors around the 
first sharp diffraction peak go , we find a strong (and inverted) response to structure at wavevectors 
below this peak: an increase in S{qo/2) lowers the degree of arrest over a wide g-range. This strong 
sensitivity to 'caged cage' packing effects, on length scales of order 2d, is much weaker in attractive 
glasses where short-range bonding dominates the steric caging effect. 

PACS numbers: 64.70.Pf, 61.20.-p 



The mode-coupling theory (MCT) of the glass transi- 
tion is a dynamical theoryof viscous slowdown and arrest 
in glass- forming liquids 0, |3| ■ It gives a good descrip- 
tion of density relaxation in molecular and (especially) 
colloidal glass-formers 3]. However, the mathematical 
structure of MCT does not readily suggest a description 
in terms of the motion of individual particles, and does 
not lend itself easily to connections to other approaches, 
e.g. free energy landscapes (see e.g. |4j). Only the cage 
effect is recognized as a major ingredient to MCT 
glassy arrest within the theory is often driven by struc- 
tural features on the nearest-neighbor lengthscale - the 
structure factor peak value S{qo). 

One major challenge is to understand how the struc- 
ture of the glass-forming liquid affects the arrested state. 
We address this issue here, focussing on the role of 
density modulations at lengthscales beyond the near- 
neighbor spacing that defines the primary 'cage'. Such 
modulations were seen in recent experiments and simula- 
tions on glass-formers that include silica melts, aromatic 
compounds, and gelation models IESSSEjOI; their 
effect on the glass is poorly understood. 

Our approach mirrors one used in superconductivity 
theory (following ^^), in which, rather than compute 
the critical temperature for a given phonon spectrum, 
one sees how Tc responds to small changes in this spec- 
trum Arrest in glasses is characterized not only 
by the temperature at which ergodicity is lost, but by 
a wavelength-dependent nonergodicity parameter F{q) 
(defined below). This describes the degree of arrest at 
different length scales, and gives the contribution to elas- 
tic scattering of the glass. We therefore choose to study 
the response of F{q) to small changes in the static struc- 
ture factor S{q'), since F{q) is well-described by MCT 
even when there is significant error in the location of the 
glass transition in the phase diagram . 

Applying this method to repulsion-dominated glasses, 
we find an unexpected sensitivity to structural features 
at wavevectors below that of the main peak in S{q). This 
is suggestive of a 'caged cage' mechanism, in which struc- 
ture just beyond the nearest neighbor shell is implicated 



in arrest. Modulation on this length scale, creating a 
positive feature in S below the peak, leads to reduced 
freezing (smaller F) over a wide g-range: such modula- 
tion thus has a cage-softening effect. The same effect is 
almost absent in glasses stabilized by short-range attrac- 
tions, reinforcing the idea that a distinct arrest mech- 
anism (bonding not caging) is at work in such glasses. 
Hence, although in MCT S{q) is the sole input to the 
glass transition, studying SF{q)/6S{q') can reveal mech- 
anistic subtleties that might otherwise be hard to discern. 

To begin, recall that MCT describes the relaxation of 
the Fourier modes (at wavevector q) of the time autocor- 
relations of the density fluctuations Pq{t) 
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In the hquid, F{q) = hmt^oo ^q{t), is zero. MCT how- 
ever predicts an idealized glass transition Q. Here, the 
relaxation time of the density fluctuations diverges and 
they never decay completely, leading to a nonzero F{q). 
The larger F{q), the stronger the arrest at wavevector q. 

The structure of the glass-forming liquid enters MCT 
through the static structure factor S{q) ~ (/Og(0)Pg(0)) 
which, apart from the mean number density n, is the only 
input needed. F{q) satisfies 



F{q)^S{q)M[F]{S{q)-F{q)), 



(2) 



which follows from the t — > oo limit of the MCT equa- 
tions of motion for correlators. Using a standard approx- 
imation for static three-point correlations, and factoring 
four-point dynamic ones, we have 

M[F] = ^J^'k V{q, kfF(k)F{p) (3a) 

V{q, k) - (eg • fc)c(fc) + (e, • p)c{p) , (3b) 

where dk = dk/{2Ti), = q/\q\ is the unit vector in 
the q direction, and p = q k. The direct correlation 
function c{q) obeys nc{q) = [l — S{q)^^] . 

MCT gives a dynamical equation for ^q{t), in which 
M[^q{t)] is a memory kernel that produces a frictional 
feedback mechanism, thought to be directly associated 
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with cage formation 0. This becomes stronger, eventu- 
ally leading to the glass transition, as static correlations 
(governed by c{q)) increase. 

We now consider only parameter values where F(q) 
and S{q) change smoothly with n and c{q). (This rules 
out the glass transition line itself; we return to this else- 
where O.) Let us calculate SF{q)/SS{q'). This tells us 
how the degree of arrest at q responds to a small change 
in S{q') at Equations © can be written in terms of 
q, k and p Considering small changes SF{q) and 

6S{q'), we obtain an equation of the form 



q'd/2n 
2 3 4 



dkA{q,k)SF{k) = J dkB{q,k)6S{k) 



(4) 



where A and B are found from ^ and JSJ. The cal- 
culation of SF{q)/SS{q') now reduces to the numerical 
inversion of A(q, k). 

All quantities G{q, k) in Q separate into singular and 
regular contributions as: 



G{q,k)^G{qUns6{q-k) + G{q,k),. 



eg ■ 



(5) 



The singular contribution {6 F{q)/SS{q')) sing = 
F{q){2S{q) - F{q))/ Siqf6{q - q') gives the same- 
mode response to a change in S{q). It ensures that 
the broad features of a change in S{q), such as a shift 
in peak positions due to compression, will reappear in 
F{q), although damped at higher q. 

In contrast, the regular part of SF{q)/6S{q') di- 
rectly embodies the mode coupling concept - a change 
in S{q') at q' affects F{q) at all q. In Fig.Q] we 
show {5F{q)/SS{q'))reg for a system of hard spheres 
just within the glass. Fig. 12 shows the corresponding 
{SF{q)/6g{r))^cg, which gives the response of F{q) to 
changes in the real-space radial distribution function 
g{r) = J qsm{qr){S{q) - 1) dg/(47r^rn). 

The dark patches of Fig.Q] correspond to 
{SF{q)/6S{q'j)rcg < 0, where a density modulation 
SS{q') > will reduce F{q), softening the glass. In hght 
areas, {6F{q)/SS{q'))rcg > 0, giving the opposite effect. 
The oscillations of the functional derivative are closely 
aligned with those of F{q) and S{q'), with period near 
to the hard-core wavenumber qg = 2Tr/d (cf. top panel 
and the constant-q cut in the bottom panel of Fig.^. 
We see that SS{q') > at q'd/27r ~ 1 will increase F{q) 
for (almost) all q. This confirms that strengthening the 
first peak in S{q') promotes arrest. Less known is that 
this also holds for q'd/2iT — N with N integer. 

On top of the oscillations, there are two remarkable 
features of SF{q)/SS{q'). First, Fig. shows a broad 
modulation envelope, maximal at relatively high qd/2'K « 
q'd/2-K K, 3. The underlying period go = 2i: / d means that 
this feature becomes localised ar r « d in the real-space 
result (Fig.[2Il: it signifies strong sensitivity to changes in 
structure within a region of width « 0. 2d just beyond the 
hard-sphere diameter. Its g-space width scales linearly 
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FIG. 1: Susceptibility of the hard-sphere-glass structure F[q) 
to changes in the liquid structure factor S{q'). The mid- 
dle panel shows the non-singular response {SF{q) / 5S{q'))ieg 
(see text); the bottom panel a cut through qd — 19 (d: 
diameter of the spheres). Top panel: S{q') used in the 
calculation (Percus-Yevick approximation at volume fraction 
95 « 0.5236). 
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FIG. 2: Real space result {5F{q)/Sg{r))reg corresponding to 
Fig.0 Upper panel: the radial distribution function g{r) used 
in the calculation. 
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FIG. 3: Constant-g cut through {SF{q) /SS{q'))icg at the max- 
imum ampUtude for a Lennard- Jones system with PY closure. 
Temperature T — e/ks, with e = 1 the well depth, and den- 
sity p = 1.085/(7'', with a — 1 the molecular radius. 
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FIG. 4: Constant-g cut through {SF{q)/SS{q'))rag, at the max- 
imum peak height for a square well system at packing fraction 
if = 0.5341 and attraction depth e = J.OkBT. The width of 
the attraction 5 — 0.03. 



with n; we associate it with the localization length of 
the system. Since SF{q)/5g{r) > here, short-range 
bonding, which increases g, also increases F. 

The second striking feature of 5F{q) /SS{q') is a deep, 
long trough at q' w qo/2 — n/d (Fig.P). This tells us 
that adding 6S{q') > at low q' will reduce arrest over a 
wide range of q. This feature gives a modulation with pe- 
riod « 2d in the real-space derivative (Fig.|2), helping to 
create an emphatic sign change at r ~ d, such that that 
sgB.{6F{q) /Sg{r)) = ±1 for r ^ d. This has two conse- 
quences: first, raising g(r) for r < d, e.g. by softening 
the core of the potential, reduces F{q) and weakens the 
glass, as expected. Second, the strongly positive region 
(at r > d) shows sensitivity to changes in g{r) beyond the 
nearest- neighbor shell. Fig. [21 shows that the strongest 
sensitivity is at r,„in ^ 1.5d, where g{r) has a minimum; 
filling in this minimum strongly favors arrest. Adding 
weight at g{rmin) could be achieved by promoting 'inter- 
stitial' particles between the first and second neighbor 
shells, strengthening the cages. This 'caged cage' effect 
emphasizes that the MCT glass transition is a collective 
phenomenon, and not just nearest-neighbor crowding, as 
many pictorial descriptions suggest. At the same time, 
the trough in Fig. ^ shows that density modulations at 
q' ~ T^/d, as might arise through incipient microphase 
separation 0|, strongly soften the glass. (Such modula- 
tions might allow more local 'rattle room' to develop.) 

These results are not specific to the hard-sphere sys- 
tem: to demonstrate this, we repeat the results of Fig. 
for a Lennard- Jones system of diameter a and attraction 
e in Fig. 121 (only a fixed-q cut is shown for simplicity) . The 
same features - a broad envelope peaked at qa/{2T:) « 3 
and a deep well at qa /{2Tr) w 1 - are found even in this 
system where attractions typical of many molecular glass 
formers are present. A clear low-q trough also remains in 
systems with rather stronger attraction, including hard- 
core Yukawa and broad 1 (i square- well systems. 

The occurrence of the trough can be rationalized by 



noting that the vertex V in Eq.O involves S{q) mainly 
through c{q). For q ^ qo, < S{q) <C 1 at high densities, 
as expected from the low compressibility of glass-forming 
liquids; thus, Sc/SS = \/{nS{q)'^) is large, and and one 
can expect 5V/5S = 5V/5c ■ Sc/SS to be large as well. 
Note that c(g) < 0, so agn{SV/Sc) = -1, and that M[F] 
in Eq.[21is a positive fimctional, leading to the expecta- 
tion SF{q)/SS{q') < for low q' . 

However, both the trough, and the consequences of it 
just outlined, are almost wiped out if strong and very 
short-ranged attractions are introduced into the system 
e.g. by supplementing the hard-sphere interaction with 
a narrow square well (Fig. 0J. As one enters the region 
of attraction-dominated glasses (see e.g. 0, ^3)' the 
trough shrinks in depth by a factor of around 100. The 
broad envelope moves to very high q, corresponding to 
a very strong localization of the particles near to con- 
tact. Fig. 13 shows the qualitative change for a constant- 
q cut through {SF{q)/Sg{r))^cg- The hard-sphere results 
for r K d (solid line) shows two main features: a peak 
of width « 0.2d corresponding to the q-space envelope 
(localization length), and broad peak of width « d con- 
nected to the trough. In the square-well results (dashed 
line) the first peak has shrunk in width to ~ (5, where 
S is the attraction range. The second peak, describing 
the response to changes within the hard core, has almost 
vanished. This reinforces the view that, in such glasses, 
arrest is primarily determined by bonding (enhancing the 
peak in g{r) close to contact (Fig. El) rather than by 
the steric confinement effects seemingly responsible for 
the trough. We emphasize that the attraction needed 
to remove this trough is exceptionally strong and short- 
ranged, and could not (to the best of our knowledge) be 
realized in molecular systems. 

Bearing this last point in mind, our study of SF/SS 
may have several wider implications. Indeed there are 
many glass-forming systems showin g a. 'pre-peak' or 

structural feature at g < go 0,0,1111113 
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FIG. 5: Constant-g cut through {5F{q) / 5g(r))^cs at the max- 
imum peak height for the square well system above (dashed 
line) . Full line shows cut through hard-sphere system of Fig. 

m 

Consider, for example, recent experimental and sim- 
ulation work on sodium silicate melts just above their 
MCT glass transition Here, the addition of sodium 

produces a pre-peak at q < ijo, and the glass is indeed 
weakened. This is driven by the formation of sodium-ion 
channels, but the general change of F{q) with increasing 
pre-peak strength is nevertheless in line with our results. 
Also, adding a smaller species to a hard-sphere system 
increases S{q) for low q and in an appropriate parameter 
region produces a glass-softening effect 19, 20], again in 
line with our Fig.^ 

The effects of a pre-peak on F{q) have also been stud- 
ied in recent experiments on the glassformer m-toluidine 
; these show the cusp predicted by MCT in the temper- 
ature dependence oiF{q) at various q. However, these X- 
ray diffraction studies do not reach high enough wavevec- 
tors to address the physics of the trough in 6F/SS. Com- 
parison might be possible with higher q measurements, 
perhaps by neutron scattering. 

A pre-peak also occurs in simulations of attractive col- 
loids, where the potential was modified to suppress bulk 
phase separation and allow access to the glass transition 
[lO|. Reassuringly, our work suggests that in attraction- 
driven glasses a pre-peak does not strongly alter F{q). 

Low-g structure, associated with arrested phase sep- 
aration, may be seen in attractive glasses (high-density 
analogs of dilute colloidal gels) Ullilillillii- How- 
ever, these features are at q <^ qo/2 and our work is 
probably not relevant to them. 

In conclusion, we have shown that the derivative 
6F{q)/5S{q') is a useful probe of the sensitivity of the 
degree of arrest at one wavevector q to changes in the 
glass-former's equilibrium structure at another wavevec- 
tor q'. In repulsion-dominated glasses, there is strong 
response to structure in S{q') at q' below the first peak. 
Although formally we are dealing with infinitesimal per- 
turbations, our results suggest that adding a 'pre-peak' 
to S{q) in this region may have a strong glass-softening 
effect (even independently of changes in the transition 



point itself). This effect, associated with packing ef- 
fects beyond the nearest-neighbor shell, becomes much 
weaker in attraction-dominated glasses. Our predictions 
are non-trivial consequences of the MCT picture of the 
glass transition. Further experimental information on 
structural arrest in systems with low-g modulations could 
therefore be valuable as a test of the overall applicability 
of the MCT picture. 
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